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SUMMARY 

Thia  clamant -contribution  modal  ia  programmed  ao  that  tha 
distribution  of  Infrarad  radiation  from  a  ship  target  can 
b«  calculated  for  any  aapact  and  any  infrared  waveband. 
Tha  present  application  ia  to  simulate  the  images  formed 
by  the  guidance  system  of  an  imaging  infrared  antlshlp 
missile. 

Tha  input  data  needed  are  the  surface  temperatures  and 
amissivltias  and  shapes  of  the  main  structures  and  the 
equivalent  black  body  temperatures  of  the  aky  and  sea 
surface.  The  shape  of  each  ship  la  represented  as  a 
number  of  'boxes'  with  each  aide  and  tha  top  of  each  box 
raprcientad  by  faces  of  elementary  cubes.  Tha  positions 
and  the  centres  of  each  face  are  calculated  from  the 
positions  of  three  rows  of  elementary  faces  measured  from 
plana  for  each  box. 

The  output  ie  an  array  of  numbers  representing  the  radiant 
intensity  from  tha  unscreened  facon  in  the  required 
direction  and  for  the  defined  wavoband.  This  array  can  ba 
converted  to  a  colour-codad  image  on  a  video  display  for 
testing  and  presentation  purposes. 
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1.  INTRODUCTION 

Ship*  nay  ba  attacked  by  aea>akimming  nlaaila*  and  guidad  bomba  which  hava 
lnfrarad  guidanca  ayatama.  Infrarad  and  Optical  Countarmaaauraa  Group  at 
Electronics  Raaaarch  Laboratory  ha*  baen  tasked  by  th«  Royal  Australian  Navy 
to  lnvaatlgata  th*  vulnarabil ity  of  ships  to  infrarad  weapons,  and  to  study 
the  af factivanaaa  of  possible  countarmaaauraa  such  as  signature  reduction  and 
the  uaa  of  decoys.  Guidanca  systems  of  th*  imaging  type  are  likely  to  ba 
especially  difficult  to  counter  by  decoys,  because  th*  spatial  distribution  of 
radiation  from  th*  target  in  th*  guidance  system  waveband  must  b*  approximated 
by  th*  dscoy,  Good  representations  of  th*  spatial  distribution  of  radiation 
from  targets  are  therefore  needed  for  countermeasures  studies. 

This  report  describes  a  target  modal  for  uaa  with  simulations  of  imaging 
infrared  ayatama  with  apacifiad  wavebands  which  may  ba  used  together  with  a 
missile  modal  and  an  atmospheric  attenuation  modal  to  study  tha  #f fsetivansss 
of  infrared  countermeasures, 

Tha  target  model  in  outlined  in  Ssotion  2,  and  described  in  detail  in 
Section  3.  Samples  of  th*  displayed  output  ere  given  in  figure  1, 


,3.  OUTLINE  OF  METHOD 

The  objective  of  th*  simulation  is  to  provide  th*  distribution  of  target 
radiant  intersity  presented  perpendicular  to  tha  eight  line  from  tha  minalla. 
Thia  ie  achieved  by  considering  smell  element*  of  th*  ship  surfaces  in  turn, 
and  calculating  their  position*  and  whether  or  not  they  ar*  visible  from  th* 
mlseile.  If  thay  era  not  screened  by  other  element*  the  position*  of  their 
projection*  on  a  plan*  perpendicular  to  tha  sight  line  from  tha  mlsaila  (th* 
'output*  plana)  ara  calculated,  and  amall  areas  centred  on  thee*  projections 
ar*  recorded  a*  parts  of  tha  projected  image.  At  the  earn*  time  these  areas 
ar*  aaalgncd  eh*  temperatures,  amiseivltlae  and  inclinations  to  th*  vertical 
associated  with  th*  elements  concerned,  which  ar*  those  assigned  to  th* 
surfaces  of  th*  major  atructuree  of  tha  ship. 

Th*  shape  of  th*  ehlp  ia  defined  in  term*  of  a  email  number  of  boxes,  with 
eaeh  face  of  each  box  uonildavad  as  made  up  of  a  large  number  of  elementary 
faces.  These  elementary  faces  are  th*  vie  ibis  faces  of  the  large  number  of 
elementary  cube*  representing  the  box  shape.  The  calculation*  for  tha 
projections  ar*  beaed  on  the  element  positions  relative  to  ship  axes,  and  th* 
position*  of  th*  misnil*  end  ship  and  the  heading  of  th*  ship  in  reference 
exes  (figure  2). 

An  a  simple  example,  figure  3  represents  a  target  considered  as  on*  box  only. 
The  shape  data  are  the  i  coordinatea  of  th*  centres  of  the  squares  in  th*  top 
row  and  bottom  row  of  th*  starboard  side;  and  th*  y  coordinates  of  th*  centres 
of  th*  aquaree  in  th*  starboard  aid*  row  of  th*  top  faoe,  ell  measured  from 
th*  first  x  position  at  tha  atarn  to  the  last  at  th*  bow.  These  three  sate  of 
data  ar*  uaad  as  described  in  Section  3.2  to  define  th*  position*  of  every 
elemental  face  in  each  face  of  th*  box. 

Th*  temperature,  emiaeivlty  and  inclination  data  ar*  the  asperate  value* 
assigned  for  the  whole  of  th*  top,  port  aid*,  starboard  aid*,  stern  and  bow. 
If  as  shown  the  box  come*  to  a  point  at  tha  bow,  there  is  no  bow  face.  Tha 
valuta  of  temperature  may  be  nominal  value*  or  values  calculated  as  described 
in  Section  3.3.  A  nominal  'ship  centre'  U  alec  defined  in  th*  ship  axle 
system,  Th*  deck*  ar*  assumed  to  ba  horliontal,  and  the  aides  vertical  unless 
otherwise  defined. 
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Thar*  r.*n  b*  several  different  *l«n*nt»  corresponding  to  any  given  point  on 
the  output,  plan*.  A  calculation  of  ranga  to  the  missile  ia  carried  out  for 
each  element  to  anaura  that  th*  on*  cloauit  to  the  missile  i*  used  to  define 
the  aurfac*  temperature  for  each  output  plana  point. 

Tha  array  of  temperatures  corresponding  to  elements  of  the  output  plane  is 
converted  to  an  array  of  radiant  intensities  by  using: 

(a)  th*  presented  area  represented  by  each  • lament i 

(b)  th*  temperature  and  emiseivity  of  each  element; 

(c)  th*  appropriate  gray-body  distribution*  convolved  with  th  required 
spectral  responses;  and 

(d)  the  inclination  of  each  element,  and  eky  and  aee  surface  equivalent 
black  body  temperatures  which  at1*  required  for  the  calculation  of  ship 
surface  radlancea. 

A  spectral  distribution  of  tha  radiant  intsnsity  is  ratainad  for  each  valua  of 
temperature  used  so  that  tha  rssults  may  bs  further  convolved  with  an 
atmospheric  transmission  model  and  any  othar  tasker  waveband  whan  required. 

The  results  may  be  presented  aa  an  array  of  surface  temperatures  or  radiant 
intensities,  or  aa  a  colour-coded  image, 


3.  DESCRIPTION  OF  MODEL 


3.1  Input  required 

Tha  input  saotion  ia  tha  first  part  of  tha  flow  diagram  of  tha  modal  shown 
in  figure  4.  After  initialising  tha  variabloa,  tha  program  prompts  tha 
user  to  supply  details  of  tha  engagement  to  be  simulated.  It  ia  naceaaary 
to  supply  tha  position  and  heading  of  the  ship,  and  the  missile  position, 
all  with  raspaot  to  a  sat  of  reference  axes;  and  to  define  the  particular 
ship  which  ia  t  •  be  eimulated  and  the  typ*  at  output  display  required.  At 
this  stage  it  is  poaalble  to  limit  the  calculations  to  ona  of  two  infrared 
wavebands  (3  to  3  p  or  t  to  14  pm), 

The  data  defining  the  ahapt  and  temperatures  of  tha  designated  vessel  are 
thon  read  from  the  appropriate  file,  which  defines  ths  small  number  of 
"boxes"  used  to  represent  tha  shape,  and  the  three  arrays  of  coordinates 
defining  the  starboard  face  of  each  box  aa  described  in  figure  3.  It  also 
includes  tha  temperatures  amisstvi.tiea  and  inclinations  to  ba  assigned  to 
each  faoe  of  each  box. 

The  temperatures  usad  era  derived  from  values  messur  ,1  in  overcast 
conditions ,  Temperature  increases  caused  by  solar  ’  .sting  may  ba 
Incorporated  by  increment ing  these  tamp* returns  a  described  in 
Section  3,3. 

3,7.  Calculation* 

3.2.1  Faoe  perimeter* 

Th*  ponitlons  of  all  th*  element*  forming  boundaries  of  facet  for  eech 
box  are  caloulatad  by  reflection  and  manipulation  of  tha  input  data  rtad 
from  fila. 


The  positions  of  the  elements  in  the  perimeter  of  tlia  port  fsca  of  each 
box  are  obtained  from  a  simple  reflection  of  the  input  starboard  face 
data  in  the  x-r  plane  which  is  the  plane  of  symmetry  of  the  ship, 

The  perimeter  of  the  top  face  of  aach  box  la  obtained  from  the  input 
data  for  the  top  of  tho  starboard  face  together  with  the  aims  data 
reflected  iri  the  x-z  plane. 

The  stern  face  ie  considered  to  be  vertical  with  boundaries  defined  by 
the  poaitions  of  tha  rear  upper  and  lower  corners  of  ths  starboard  and 
port  facai.  If  tha  stern  upper  elements  of  the  port  and  starboard  aides 
ars  adjacent,  there  is  no  ..tern  face. 

Tha  perimeter  of  t.he  front,  face  is  similarly  defined  from  tha  positions 
of  tha  bow  upper  end  lower  elements  of  the  port  and  starboard  facta. 

3.2.2  Calculation  of  computing  loop  limits 

Tha  computations  on  tha  surface  el  ament r  to  determine  their  projected 
positions  on  ths  cutout  plane,  and  thair  radiance  values  are  carriod  out 
in  "naated"  iterations.  Tha  outer  neat  coven  the  box  number,  and  tho 
next  the  faoo  number  (eg  Starboerd-1,  Port-2,  Top-3,  Stern-4,  Bow-5), 
The  next  neat  covers  ths  sections  along  tha  ship  axis  (ths  x  coordinate 
changes)  from  ths  st.ern  to  the  bow.  Ths  final  loop  is  either : 

(s)  over  the  range  of  valuta  of  height  (s)  at  constant  distance  from 
tha  x-s  plana  (y)  for  tha  vertical  tides!  or 

(b)  over  the  range  of  values  of  y  from  ont  aids  to  ths  other  of  the 
flat  top  at  constant  t)  or 

(o)  over  th#  rang*  of  valuta  of  both  y  and  s  for  tho  stern  and  bow. 

It  should  be  notad  that  the  cross  sections  of  all  (acta  porpandiculsr  to 
the  a  axis  are  represented  as  althsr  vertical  or  horisontal  lines.  This 
limitation  doss  not  affact  the  appearance  of  tha  images  to  any 

noticeable  degree.  However  the  inclinations  of  the  faoes  to  the 

vertical  era  taken  into  account  for  calculating  tha  contributions  of 
reflection  from  sky  and  sea  surface  to  total  radiant  intensities  (eoo 
Appendix  II). 

3.2.3  Element  mapping 

The  rotation  constant*,  which  depend  on  the  heading  of  the  ship  and  its 
poaltlon  relative  to  the  missile,  and  which  are  independent  of  position 
on  the  ship,  ars  calculatsd  in  a  subroutine  outside  the  iteration  loops 
in  which  the  mapping  of  element  positions  onto  thu  output  plane  are 
performed, 

The  mapping  ie  carried  out  by  converting  thr  x,y,s  coordinates  of 
Individual  element*  int.o  three  coordinates  ,  y^,  z^)  in  a  system 

where  X_  is  ths  direction  from  observer  to  shi,'  Y_  has  ths  direction  of 
P  P 

v  x  Xp  where  v  is  the  upward  vertical,  and  Z^  complete*  th*  right  handed 

orthogonal  sat.  2  Y  is  the  output  plane  and  y  ,  s  the  coordinates  of 
P  P  F  P 

the  projected  point  in  this  plane.  The  equation  for  thia  conversion  ie: 
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where  A i ,  Aj  are  tho  familiar  rotation  matrices  (hoc  Appendix  I)  and  if< 

and  6  are  rotations  which  are  derived  in  Appendix  1, 

As  oach  element  is  mapped  onto  the  output  plane,  x^  provides  a  measure 

of  the  range  from  the  element  to  the  missile.  The  temperature  assigned 
to  esch  position  in  the  output  plane  is  that  associated  with  tho  element 
cloaest  to  the  missile  which  has  bean  mapped  onto  that  position.  By 
this  means  the  masking  of  any  part  of  the  ship  by  any  other  part  1b 
incorporated  in  the  model. 

Owing  to  rounding  errors  in  the  computations  it  was  found  that  small 
sections  of  tho  rotated  imago  were  blank  or  showed  sections  which  should 
have  buen  masked.  This  was  overcome  by  assigning  the  temperaturn  of  any 
element  of  the  output  plane  to  the  next  in  increasing  value  of  and  b 

(and  a  fourth  element  to  complete  the  square).  Those  additional 
element*  were  normally  overwritten  by  subaequont  mapping,  but  the  effect 
waa  to  fill  in  the  minor  gape,  A  further  effect  of  this  technique  was 
to  incronso  tho  sisse  of  each  mapped  face  by  one  olemeni;  in  the 

direction  and  one  element  in  tho  direction.  If  tho  element  sine  is 

small  compared  with  ths  fees  slr.e  tills  affect  is  of  no  consequunco.  The 
effect  can  however  be  eliminated  by  reducing  the  input  data  values  by 
one  unit  in  each  direction, 

3.2.4  Calculation  of  radiatod  and  reflected  radiant  intensities 

Ths  radiant  intanalty  of  each  element  of  the  output  piano  array  depends 
on  the  presented  area  A  it  represents  (which  Is  the  same  for  all 

elements  in  ths  output  plans),  the  temperature  (T°K)  and  smlaaivity  (t) 

of  the  whip  olement  it  represents,  and  the  infrared  waveband  chosen. 
The  radiant  intansitios  (W/nt)  are  calculated  using  the  Planck  formula: 

Xj  . 

,  C,  A  *  B.  dX 

HX»-X,  "Vj  "xpCCi/VOM  (z) 

X  i 

whore  C*  and  C|  are  the  first  and  second  radiation  constants 

3,74  x  lo’11  W  cm',  and  1.438  cm8K  respectively,  \  is  wavelength  in 

centimetres,  and  is  the  normalised  spectral  sensitivity 

characteristic  of  the  sensor  from  wavelength  X,  to  A(. 

In  addition  to  the  enurgy  radiated  from  tho  surface,  a  proportion  of  the 
energy  originating  in  tho  hemisphere  it  subtends  is  ref  looted ,  Tho 
additional  radiant  intensity  from  this  soure.n  is  given  byi 
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A  (1  -  t)  Et 
JREFX,-X,  "  * 

An  expression  for  F,^,,  the  energy  received  from  the  background  by  the 

element,  In  derived  In  Appendix  II.  The  cnee  coimidersd  Is  for  a 
uniform  sea  background,  a  sky  background  radiance  which  is  a  function  of 
elevation  angle  only,  and  a  perfectly  diffuse  surface  of  the  element  A. 

In  order  to  calculate  the  lrradiance  at  the  misailo  from  each  element  it 
is  also  nacessary  to  incorporate  the  atmospheric  transmission  (t^)  in 

the  integration  of  equation  (3).  At  present  this  is  pot  formed  by  using 
Lowtran(ref .  1) . 

Since  there  are  in  general  only  about  ten  boxes  each  with  five  faces, 
there  are  only  about  fifty  temperatures  to  consider.  The  radiant 
intensity  and  irradianca  calculation!  ora  performed  outside  the 
iteration  loops  and  only  once  for  each  value  of  temperature  assigned  in 
the  model  to  avoid  repeating  the  calculations  for  each  element. 

3,3  Output  presentation 

The  output  may  bo  presented  as  a  printed  array  of  the  values  of  surface 
temperature,  radiant  intensity,  or  irradlanr.ej  or  as  an  array  of  symbols 
with  a  key  for  conversion  to  thu  necessary  values, 

A  brack  and  white  shadod  image  can  be  used  to  represent  the  spatial 
distribution  of  tempnraturo,  radiant  intsnslty  or  irradlunce  with  an 
appropriate  key,  and  colour  can  ba  used  Instead  of  or  as  well  as  shading. 

The  type  of  presentation  required  is  selected  by  the  user  at  the  start  of 
tha  simulation. 


4.  PRESENT  APPLICATION  OF  THE  MODEL 

The  target  model  is  intended  to  be  part  of  a  simulation  of  mlssilo  attacks  on 
ships,  during  which  infrared  countermeasures  may  be  used.  Sr-  far  the  goometry 
of  nn  FFG  Frigate  (HMA8  Canberra),  a  Destroyer  Escort  (HMAS  Torrsns)  and  a 
Guided  Missile  Destroyer  (HMAS  Perth),  and  a  Coastal  Minesweeper  (Mill)  have 
been  modelled,  Surface  temporaturas  have  bean  derived  from  trials  in  which 
thermovision  cameras  recorded  the  infrared  images  of  some  of  these  ships. 

The  background  is  represented  as  a  constant  lrradiance  from  thu  sea  surface, 
and  an  lrradiance  equivalent  to  that  from  a  black  body  nt  ambient  temperature 
(with  no  atmospheric  attenuation)  for  the  horlxon  sky.  The  lrradiance  from 
the  sky  reduces  to  noro  linearly  with  elevation  angle  up  to  90®,  This 
represents!  ion  requires  rnflnement.  for  representations  of  specific 
environments,  The  position  of  ths  waterline  cm  the  ship  rulative  to  the 
horlaon  is  calculated  from  the  height  of  the  missilo,  tha  range  of  the  ship, 
and  the  radiUN  of  the  earth  (Appendix  111), 


S ,  MODEL  DEVELOPMENT 
3,1  Irradianca  from  exhaust  plum os 

For  missiles  using  the  infrarad  waveband  which  covers  thu  carbon  dioxide 
emteniona  near  4,4  pm,  the  engine  exhaust  plumeu  can  be  significant 
sources,  This  is  particular  ly  so  for  tho  FFC)  vessels  which  have  gas 
turbine  engines. 
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A  series  of  trial*  were  carried  out  for  various  purposes  on  an  FFG  in 
March  1963(ref.2),  The  video  records  of  plums  signatures  are  being 
analysed  in  an  attempt  to  derive  empirical  models  of: 

(a)  irradiance  from  the  plumes  as  a  function  of  range,  aspect,  and 
engine  condition!  and 

(b)  plume  shape  aa  a  function  of  wind  speed  and  direction,  ship  speed 
and  direction,  aspect,  and  angina  condition. 

During  the  asms  triala  period  limited  measurements  were  also  made  on  a 
daatroyar  escort  with  a  steam  turbine. 

Examples  of  8  to  14  and  3  to  5  dm  images  of  HMAR  Canberra  are  shown  in 
figure  5.  The  plume  ia  the  dominant  feature  of  the  3  to  5  pm  signature  of 
the  FFQ  at  the  range  of  the  measurement,  The  images  shown  wars  produced 
with  the  analysis  system  described  in  reference  .1. 

It  ie  hoped  that  it  will  ba  possible  to  incorporate  a  three-dimensional 
representation  of  FFQ  and  DF.  plumes  in  the  target  model  following  analysis 
of  trials  data. 

3.2  Background  representation 
(1)  Sea  background 

The  radiance  of  the  sea  surface  is  partly  grey-body  radiation  at  the  see 
temperature  and  partly  reflection  of  the  sky,  Since  sky  rsdisnes  is  in 
gsnsrsl  a  function  of  elevation  angle,  the  aoa  radiance  will  be  ii 
function  of  the  declination  of  the  eight  line.  Sea.  roughnose  may  bo 
considered  as  the  variation  of  slops  of  areas  of  surface  with  certain 
spatial  and  temporal  corveUtion  characteristics.  Variation  of  aurfacu 
angle  will  produce  variations  of  the  reflected  component  of  sea 
radiance,  To  produce  a  target  model  which  incorporates  sea  clutter,  the 
surface  shape  variation*  must  be  modelled  so  that  the  spatial  and 
temporal  changes  are  derived  as  functions  of  local  wind  speed  and 
direction  (and  possibly  other  parameters  sf footing  sea  surface  shape), 
Such  models  do  exist,  for  oxample(rcf ,2,3),  and  considerable  work  has 
boon  carried  out  to  estimate  the  accuracy  of  modol  predictions , 

Fluctuations  of  background  signals  are  often  the  median  lams  for  limiting 
tile  detectability  of  targets  with  particular  systems,  and  may  bu  even 
more  important  in  a  countermeasures  environment.  It.  is  recommended  that 
a  dynumlc  sun  background  modol  bo  dovelopod  for  use  with  the  present 
target  model. 

53  Derivation  of  ship  surface  tempo ratures 

5.3,1  Equilibrium  temperatures  of  solar  heated  surfaces 

To  obtain  initial  estimates  of  temperatures  of  the  sunlit  portions  of 
ship  surfaces,  ths  plates  may  be  considered  as  sufficiently  thin  or 
thermally  conducting  to  allow  both  eurfacea  to  be  at  the  same 
temperature,  and  to  have  negligible  thermal  mass.  The  plates  are 
irradiated  by  the  sun  on  one  side,  and  exchange  host  by  radiation  and 
convection  with  the  aurroundings  on  both  sides, 

Ths  thermal  input  powers  per  unit  area  of  plate  surface  ere: 
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(1)  Ficoiu 

(2)  spoT,‘ 


Th*  solar  irrsdlanco  on  thn  place  where  F  1h  the 

s 

sea-levul  normal  solar  irradiance,  and  u  is  the  angle 
of  the  normal  of  the  plate  to  the  solar  sight-lino. 

The  radiation  absorbed  from  ambient  at  absolute 

temperature  Tj  onto  aide  1  of  the  plat*  whnre  sp  1* 

the  emisalvity  of  thn  plate  surface,  and  a  is 
Stefan'a  Conatant, 


(3)  tpOTi*1  As  (2)  but  for  side  2  of  the  plate. 

Th*  thermal  output  powers  per  unit  area  of  plat*  surface  arei 

(4)  aF^coiw  Th*  aolar  radiation  reflected  from  thn  surface,  wnere 

a  it  th*  plate  surface  reflectivity  for  the  solar  sea 
lavel  spectrum. 


(5) 

yV 

The  radiation  from  the  plate 

T  into  ambient  1. 

surface 

at  temperature 

(6) 

yV 

Aa  (5)  but  into  ambient  2. 

m 

(hi+BiVi 

)(T  ‘TO  The  convection  from 

surface 

1  of  the  plat* 

into  ambient  1  whore  ht  and  Bi  are  free  and 
forced  convection  coefficients  end  Vt  is 
the  wind  velocity, 


(8)  Ou+BiViHT^Ti)  The  convection  from  surface  2  of  the  plate 

into  ambiant  2. 


Th*  thermal  power  aquation  far  equilibrium  is  therefore i 


F^cosw  +  tpoT,‘  +  ipdTi* 


sF^cosw  +  2*  oT#4  +  (h,+#,V,)(TM-T,) 


+  (li.+Bj  V, )  (T^  -T| ) 


This  equation  can  be  solved  to  provide  the  values  of  T  by  providing 
known  values  of  F(,  w,  »pl  o,  a,  T i ,  Tt ,  V,  and  V(  and  choosing 
appropriate ,  values  of  hi,  lit,  Bi  and  fit,  the  vnluas  for  the  outer 

•urface  depending  on  windnpetd  and  aspect,  It  is  expected  that 
attributing  values  of  Vj  for  net  wind  spued  across  the  various  faceH  of 

th*  ehip  will  be  difficult,  owing  to  the  degree  of  nxposur*  to  or 
ahelter  from  the  relative  wind  direction. 

It  ehould  be  noted  that  T,  ie  an  offer, tive  temporature  for  the 

hemiaphere  aubten.ied  by  a  surface.  It  can  bt  obtained  by  measurement 
using  a  radiometer  with  sufficiently  wide  waveband  to  covet  all 
significant  energy  inputu.  Alternatively  Tj  could  be  calculated  front 

atmospheric  data  to  Include  tingle  scattering  and  radiation  from  the 
molecules  of  the  atmosphere. 
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5.3.2  Derivation  of  non-equilibrium  surfaco  temperatures 

Because  the  plates  of  the  ship  do  have  thermal  mess,  the  equilibrium 
temperatures  are  only  reached  after  some  time  of  continuous  irradiation 
under  constant  conditions.  In  general  solar  irradiation  is  not  constant 
owing  to  changes  of  sun  angle,  cloud  cover,  and  ship  heading. 

The  equation  of  Section  5.3.1  can  be  modified  to  incorporate  changes  of 
Fs  and  u  and  to  include  thermal  time  constants  of  plates  based  on  their 
thickness,  specific  gravity  and  specific  heat.  If  the  plates  are  thin 
and  of  metal,  thermal  conductivity  should  be  sufficiently  high  to  keep 
both  surfaces  at  vary  similar  temperatures. 

In  the  absence  of  solar  irradiation  the  outer  surfaces  of  the  ship  will 
reach  equilibrium  temperatures  between  the  inside  temperatures  (which 
include  those  of  funnels  and  engine  rooms)  and  the  outside  ambients. 

5.3.2  Development  of  temperature  prediction 

The  model  may  be  developed  by  programming  the  equation  of  Section  5.3.1 
and  choosing  the  appropriate  heat  transfer  coefficients  to  provide 
equilibrium  temperatures.  A  solar  position  model  ia  already  available 
to  provide  w  from  latitude,  longitude,  Zulu  time,  and  ship  geometry  and 
heading. 

Further  development  to  derive  non -equilibrium  temperatures  along  the 
lines  of  reference  4,  requires  research  on  plate  materials  and 
thicknesses  and  ship  internal  temperatures.  (For  the  Frigates  with 
superstructure  rf  i  inch  aluminium  the  thermal  capacity  .s 

4.1  Uhl1  °c"1.  The  corresponding  value  for  the  steel  hull  is  10.4  Dr 
8.65  depnnding  on  position,) 

Such  developments  would  allow  prediction  of  ship  signatures  after 
planned  voyages  and  could  be  valuable  in  special  military  circumstances. 

5.4  Extension  to  other  ships  and  military  targets 

It  is  intended  that  one  of  each  of  the  major  ship  types  in  the  Royal 
Australia, i  Navy  will  be  modelled  when  suitable  plans  become  available. 

Although  the  intention  ia  to  develop  methods  of  reducing  the  vulnerability 
of  RAN  ships,  it  would  also  be  reasonable  to  model  foreign  ships  to 
investigate  their  vulnerability  in  the  infrared  wavebands. 

The  method  car.  also  be  applied  to  generating  infrared  images  of  other 
military  vehicles  such  as  tanks. 


6.  CONCLUSIONS  AND  RECOMMENDATIONS 

A  digital  model  for  generating  the  infrared  images  of  ships  has  been 
developed.  Shape  data  for  four  ships  have  been  used  with  the  model  to  produce 
images  and  arrays  of  irradiances  suitable  for  inputs  to  a  missile  model. 

An  engine  plume  model  is  needed  to  complete  the  images  in  the  3  to  5  pm 
waveband.  An  improved  representation  of  background  radiances  is  required  for 
simulation  of  specific  situations.  A  detailed  model  of  the  thermal  effects  at 
the  ship  surfaces  is  necessary  to  extend  the  model  bo  that  surface 
temperatures  are  calculated  from  the  thermal  history. 
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NOTATION 


C, 

Eir 


JRKF 


"Rl 

R 


Rj 

R> 


T 


T 


Rl 


T 


s 


V 


X  1 


I! 

« 


area  (m* ) 

rotation  matricea 

normalised  apaotral  aanaltlvlty  at  wavelength  X 
flrat  radiation  constant  (W  cm*) 
second  radiation  constant  (on  °K) 
lrradianca  from  background  (V  on'1) 
irradianca  from  upper  hemisphere  on  area  A  (W) 
lrradianca  from  lower  hemisphere  on  area  A  (W) 
sea  level  normal  solar  lrradianca  (W  cm  * ) 
radiant  intensity  (W  sr"1) 
reflected  intensity  (W  nr"1) 
missile  position  in  reference  system 
range  (km) 

radius  of  Earth  (km) 
axes  subscript  for  reference  system 
axes  subscript  for  sight-line  system 
temperature  (®K) 

target  position  in  reference  system 
temperature  of  sea  surface  (°K) 
velocity  (cm  s  *) 

axes  systems 

ref lsctivity 


elemental  areas  (cm*) 


frss  convection  coefficient  CW  a"*  •jc'1) 
bow  coordinate  of  ship  element 
port  coordinate  of  ship  eleaent 
height  coordinate  of  ship  eleaent 

ship  eleaent  coordinates  in  sight-line  system 

asimuth  angle  limit  (•) 

inclination  to  vertical  (*) 

forced  convection  coefficient  (W  in'*  ®x  a) 

angle  of  reflection  (•) 

declination  of  sight-line  (•) 

declination  of  horlaon  (*) 

angular  separation  of  horiaon  and  sight -line  (•) 

amiss ivity 

Euler  rotation  (®) 

wavelength  (pm) 

aeimuth  angle  (®) 

Stafan'e  conetant  (V  sT*  °k'®) 
elavation  angle  (*) 

Euler  rotation  (®) 

angle  of  plate  normal  to  eight-line  (•) 
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APPENDIX  I 

DERIVATION  OF  ROTATION  EQUATION 
1.1  Sight-line  exee  (R2) 

All  point*  on  the  target  are  mapped  onto  i  plane  perpendicular  to  tho  sight 
line  from  the  missile  to  the  target  (thi  output  plane),  flight-line  axes 
are  in  the  direction  of  the  sight  line,  YR2  oerpendiculsr  to  the  upward 


vertical  and  to  XR2,  and 

ZR2  to  complete  ths  righe-hsndsd  orthogonal  set 

let 

XR2  ' 

<TR1  *  MR1^ITR1  ■  MR1* 

YR2  " 

(ZR1  X  XR2}/*ZR1  X  XRal 

ZR2  " 

^2  x  YR2 

where  TR1  (x^  y^i  Sj.)  and  MRJ  (x^,  yH,  *M)  are  the  position  vector*  of  the 
target  and  missile  referred  to  ths  origin  of  reference  axes  XRJ,  YR1,  2R] 
with  ZR1  ohosen  to  be  the  upward  vortical.  The  ship  is  assumed  to  be 
moving  perpendicular  to  the  ZR2  axis  with  *T  ■  0  and  with  no  roll,  pitch, 
or  yaw. 

1.2  Rotation  transform 

Every  point  x,  y,  s  referred  to  target  ax**  (X^  bow,  Y^  port,  ZT  vertically 
up)  is  expressed  in  coordinate*  (xpl  y^ ,  ap)  referred  to  sight* line  axes, 

the  Y-Z  plan*  of  which  is  th*  output  plan*  so  that  screening  of  any  point 
by  another  oan  b*  checked,  and  to  provide  tho  data  for  presentation  of  an 
image  of  th*  ship  as  ssen  from  th*  missile, 

Two  Euler  rotations  end  8  (figure  1.1)  sr*  used  to  rotate  ths  target  axes 
(T)  to  align  them  with  th*  reference  axe*  (Rl).  9  is  u  rotation  about  Z^ 

to  bring  Y^  to  Y'  ■  YR2,  perpandiuoUr  to  XR2,  and  X^.  to  X'  let 


(and  is  the  angle  between  the  alght  line  and  th*  ship  bearing  in  th* 
horizontal  plans).  0  la  a  rotation  about  tha  now  Y*  axis  to  bring  tha  naw 
X1  axis  into  lino  with  th*  XR2  exist 


(end  8  is  th*  angle  between  th*  tight  line  and  ths  horisontsl  plans).  Ths 
rotation  transform  iat 
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APPENDIX  II 

REFUND  -radiant”  intensity  from I  an I  inclined  plane  element  originating  from 

oflCKwOl/ND 

II.  1  Upper  hemisphere 

'sjvs,  ts=s.  s.r  j 

ths  area  of  the  insids  aurfac.  of  a  aph.re  of  large  radius  R,  between  two 

cepe,  which  is  subtended  by  the  plane  element  dA,  inclined  at  o  to  the 
vertical  (figure  II. 1), 

*”* '  “1"lh  <“*-'•  «i» 

dAt  dAt  (coed  coep  cosa  -  alnd  sing)  f(d) 


where:  dA,  -  R,*  coed  dd  dp  cm* 


.  1  C,  X*«  bx  dX 

*  j  exp  C,/XT(d)  w  ®r  1  °®  * 
/  X, 


The  energy  falling  onto  dA,  from  the  whole  horizontal  section  of  the  sphere 
X*  dEu»  given  by: 


dEu  -  2dA,  coed  f(d)  dd  cosd  coso  [ 


cosp  dp  - 


aind  sins  /  dp  V 


0 
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and  tha  energy  falling  onto  dAi  from  the  whole  upper  hnmliphere  ii 


Eu  ■  2dA, 


it-a 


r 

cos’  <t  sinl  f(d)  d# 

. 


n-a  _ 

* 

sina  J  cosd  *in$  l  f(t)  dd 


where  £  ■  arcoa  (tand  tuna) 

II. 2  Lower  hemisphere 

The  apparent  background  aea  radiance  ie  aaaumed  to  be  coni t ant,  An  element 
of  lea  background  dA,  la  represented  by  part  of  the  area  of  the  inalde 

surface  of  a  segment  of  a  sphere  of  large  radius  Rs  (figure  II. 2). 

If  the  element  dA,  is  at  d  declination  and  p  arimuth  (figure  II. 2),  the 
radiation  falling  onto  dAk  from  dA,  is: 

dAi(co»fl  cosp  CQ»q  -  slnd  cosp  alna)  dA,F 

R,‘ 


where 


dA,  ■  Rj*  cosp  dd  dp  cm* 
and 

.  r  C,  X  *  BXdX 
F  "  i  j  exp  c,/vr§  w  ,r  ' 

where  is  the  absolute  temperature  of  the  sea  surface. 
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Th«  anargy  falling  onto  dAi  from  tha  who  la  aagmant  at  i  ia 


L  ■  dA|  coaa  j  coad  coa'p  d$dp  -  sinct  j  Bind  coa'p  dddp  W 


and  from  tha  whola  lowar  hamisphara: 


n+2 

2  2 

7  L 


»+2  +* 
2  2 


cost  coa'p  dddp  •  a  inn 


sing  coa'p  dddp  W 


II. 3  Reflactad  intanalty 

Tha  total  anargy  in  tha  dafinad  wavabar.d  falling  onto  dAt  which  originataa 
from  background  iai 

et  -  eu  +  elv 

A  proportion  of  E^,  (1  *  c)  la  raflactad  from  tha  aurfaca,  and  if  thia 

raflaetlon  la  aaaumad  to  ba  lambartian  (diffuaa)  tha  raflactad  intanalty  at 
any  angla  Y  to  tha  normal  of  dA|  ia; 


REF  X,-X,  ■ 


dA,  eo»Y(l-t)ET 
- T  w 


Figure  II. 2  Geometry  of  lea  background  irrariiance  calculation 
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\tl  SUMMARY  OR  ABSTRACT: 

1*^  (Ifthia  It Hcurlty  ajaialflad. thi announasmtm  of  thli  upon  will  b« aimllrly  claaaiflad)  _ 

'V*  Thia  element-contribution  modal  ia  programmad  ao  that  tha  distribution 
of  infrarad  radiation  from  a  ship  target  can  ba  calculated  for  any 
aapact  and  any  infrared  waveband.  Th*  present  application  ii  to 
simulate  the  images  formed  by  the  guidance  ayatem  of  an  imaging  infrared 
antlship  misail*. 

The  input  data  needed  are  the  surface  temperatures  and  emmlaiivitias  and 
shapes  of  the  main  structure*  and  the  equivalent  black  body  temperature) 
of  th*  aky  and  sea  surface.  Th*  shape  of  each  ship  ia  represented  aa  a 
number  of  'boxes'  with  each  aid*  and  th*  top  of  each  box  represented  by 
faces  of  elementary  cubes.  Th*  positions  and  th*  cantrea  of  such  face 
ar*  calculated  from  th*  position*  of  three  rows  of  elementary  face* 
measured  from  plana  for  each  box. 

Th*  output  is  an  array  of  numbers  representing  the  radiant  intensity 
from  th*  unscreened  faces  in  the  required  direction  and  for  th*  defined 
waveband.  This  array  can  b*  converted  to  a  colour-coded  image  on  a  video 
display  for  tasting  and  presentation  purpose), 
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